Orbital eigenchannel analysis for ab-initio quantum transport calculations 
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We show how to extract the orbital contribution to the transport eigenchannels from a first- 
principles quantum transport calculation in a nanoscopic conductor. This is achieved by calculat- 
ing and diagonalizing the first-principles transmission matrix reduced to selected scattering cross- 
sections. As an example, the orbital nature of the eigenchannels in the case of Ni nanocontacts is 
explored, stressing the difficulties inherent to the use of non-orthogonal basis sets 
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I. INTRODUCTION 

Ab initio quantum transport calculations for 
nanoscopic conductors like molecular junctions and 
nanocontacts have become standard^, and replace, to 
a large extent, earlier phcnomcnological or parametrized 
procedures for quantum transporls^. A widely ac- 
cepted starting point is the Landauer approach which 
describes the electron transport as an elastic (and thus 
phase-coherent) scattering process of non-interacting 
quasi-particles. In this approach the conductance of 
a nanoscopic region is determined by the quantum 
mechanical transmission probabilities of the incoming 
transport channels 5 . The transmission matrix is calcu- 
lated by mode-matching of the incoming waves in one 
lead with the outgoing waves on the other lead via the 
intermediate scattering region 1 or, equivalently 6 , from 
the one-body Green's function of the scattering region 
connected to the metallic leads^. The latter is known as 
the non-equilibrium Green's function (NEGF) approach 
since a finite bias voltage can also be taken into account 
when calculating the electronic structure of the junction. 
One of the advantages of the NEGF approach is that 
it does not need to be implemented from scratch since 
one can make use of standard ab initio codes for the 
computation of the electronic structure. 

It is often useful to decompose the total conductance 
into the contributions of the transport eigenchannels, 
first introduced by Biittiker 7 .. These are defined as the 
linear combinations of the incoming modes in a lead that 
do not mix upon reflection on the scattering region and 
present a unique transmission valued. The decomposi- 
tion of the measurable total transmission in terms of 
the transmissions of these eigenchannels simplifies con- 
siderably the interpretation of the results. Knowledge 
and analysis of the eigenchannel wavefunctions would, in 
turn, allow one to make predictions regarding the behav- 
ior of the conductance upon distortions of the geometry 
or other perturbations of the scattering region^. Unfor- 
tunately, in the NEGF approach it is not straightforward 
to extract the orbital composition of the transport eigen- 
channels. Only the eigenchannel transmissions can be 
obtained easily in the NEGF approach from the non- 
negligible eigenvalues of the transmission matrix. How- 
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FIG. 1: Sketch of the scattering problem. L: Left lead. D: 
device. R: right Lead. 0: cross-section of interest. 



ever, the associated eigenvectors turn out to be useless as 
obtained. The reason is that these eigenvectors contain 
the contributions to the eigenchannel wavefunctions of 
the atomic orbitals at one of the borders of the scatter- 
ing region immediately connected to the leads, but not 
inside. 

Here we present a method for analyzing the orbital con- 
tributions to the transport eigenchannels at an arbitrary 
cross-section of a nanoscopic conductor by calculating 
the transmission matrix projected to that cross-section. 
Our approach generalizes previous work by Cuevas et ali^ 
for tight-binding-type Hamiltonians to non-orthogonal 
atomic orbitals basis sets as those commonly used in 
quantum chemistry packages. An alternative approach to 
investigate the contributions of certain atomic or molec- 
ular orbitals to the conductance consists in directly re- 
moving the respective orbitals from the basis seljiS. 



II. METHOD 

First we divide the system under study into 3 parts: 
The left lead (L), the right lead (R), and the intermediate 
region called device (D) from now on, where only elastic 
scattering takes place. Figure^ shows a sketch of a nano- 
constriction connecting two bulk leads. We assume that 
the leads are coupled only to the constriction but not to 
each other. The Hamiltonian describing this situation is 
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then given by the matrix 
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Since many density functional theory (DFT) codes work 
in non-orthogonal basis sets, we also allow explicitly for 
overlap between atomic-orbitals given by the following 
overlap matrix: 
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The standard approach to calculate the conductance 
is to calculate the self-energies of the leads from the 
Green's functions of the isolated leads, i.e., for the left 
leadS L (S) = (H DL -£SDL)gL(£)(H L D-£S LD ) where 
gh(E) = (ESl — Hl) _1 is the Green's function of the iso- 
lated left lead and analogously for the right lead. From 
this we can calculate the Green's function of the device: 

G D (E) = (ES D - H D - S L (E) - V R (E))-\ (3) 

which, in turn, allows us to calculate the (hermitian) 
transmission matrix 

T(E) = r^Ey^G^T^G^E^iE) 1 / 2 , (4) 

where T L = i(E L - £[) and T R = i(S R - S^). Typ- 
ically the leads are only connected to the left and right 
borders of the device and are sufficiently far away from 
the scattering region so that they can be described by a 
bulk electronic structure. From the structure of eq. (QJ 
it follows that only the submatrix of T representing the 
subspacc of the device immediately connected to one of 
the leads are non-zero. Thus the eigenvectors obtained 
by diagonalizing T only contain the atomic orbital con- 
tributions to the eigenchannels at one border of the de- 
vice region but not at the center where the resistance is 
ultimately determined. 

To investigate the orbital nature of the eigenchannels 
at an arbitrary part of the device we can simply calcu- 
late the transmission matrix associated to this part. By 
choosing this region to be a cross-section, like that indi- 
cated in Fig. ^ current conservation guarantees that the 
so-calculated conductance is approximately equal to the 
conductance calculated from the transmission matrix of 
the whole device. We want to emphasize here that this is 
really only approximately true for a Hamiltonian beyond 
the tight-binding approximation since hoppings between 
atoms on both sides beyond the selected region are ne- 
glected. Of course this approximation becomes better 
the thicker the chosen cross-section is. We proceed by 
further subdividing the device region. The cross-section 
of interest will be referred to as while the regions on 
either side will be denoted as 1 and r, respectively: 

(hi hm hi r \ 
hoi h h 0r (5) 
h r] h r0 h r J 



(Sl Sio Sl r \ 
soi s s 0r . (6) 
S r l S r0 S r / 

As mentioned above we will neglect the hoppings (and 
overlaps) between the left and right layers outside the 
region of interest so we set hi r = h r i = sir = s r i = 0. 
With this approximation the Green's function matrix of 
the cross-section can be written as 

G (E) = (Es - h - - K(E))- 1 - (7) 

The self-energy matrices representing the coupling to the 
left and right lead, = (h ol ~ Es ol )g[(E)(h w ~ Es l0 ) 

and V r (E) = (h 0r - Es 0l )g' r (E)(h r0 - Es t0 ), are given 
by the Green's function of the left layer 1 connected only 
to the left lead L and the right layer r connected only to 
the right lead R, respectively: 

g[(E) = (Esi-hi-^iE))- 1 (8) 
g' T (E) = (E Sl .-h 1 .-S R (E))- 1 . (9) 

The reduced transmission matrix (RTM) with respect to 
the chosen cross-section is now given by 

t'(e) = r;(s) 1 / 2 G (^)r;.(^)G (s)r((^) 1 / 2 (10) 

with r; = *(E{ - £( f ) and r; =*(££- S; 1 ). Diago- 
nalizing T(E) now yields the contribution of the atomic 
orbitals within the cross-section to the eigenchannels. 

III. RESULTS AND DISCUSSION 

In the following we apply the above described method 
to analyze the orbital nature of the conducting channels 
of Ni nanocontacts which have recently attracted a lot of 
interest because of their apparently high magnetoresis- 
tive properties^. We consider the nanocontact to con- 
sist of two ideal pyramids facing each other along the 
(001) direction and with the two tip atoms being 2.6 A 
apart. Bulk atomic distances (2.49 A) and perfect crys- 
talline order are assumed for each pyramid. Just as in our 
previous work on Ni nanocontacts^ we perform ab initio 
quantum transport calculations for this idealized geome- 
try. To this end we use our code ALACANT (ALicante 
Ab initio Computation Applied to NanoTransport). The 
electronic structure is computed at the DFT local spin 
density approximation level with a minimal basis set and 
the electrodes are described by means of a semi-empirical 
tight-binding Bethe lattice model. 

As indicated in Fig. Hwe calculate the RTM T'(E) for 
one of the tip atoms of the contact (labelled with 0) and 
diagonalize it to obtain the eigenchannels and the cor- 
responding transmissions projected on the tip atom. In 
Fig. |3 we compare the individual channel transmissions 
calculated on the one hand from the full transmission ma- 
trix (FTM) T(E) and on the other hand from the RTM 
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FIG. 2: Sketch of Ni nanocontact consisting of two pyramids 
facing each other along the (001) direction with the two tip 
atoms forming a dimer bridge. The device region (grey cir- 
cles) consists of 28 Ni atoms and the left and right electrodes 
(empty circles) are modelled by Bethe lattices with appropri- 
ate tight-binding parameters to reproduce Ni Bulk DOS. 



(a) majority chan. 



(b) minority chan. 1 (c) minority chan. 2,3 




FIG. 3: Transmission functions of open transport channels 
for the Ni nanocontact sketched in Fig. [5] as calculated from 
the FTM T(E) (solid line) and from the RTM T'(E) (dashed 
lines), (a) shows the only contributing M channel and (b)-(c) 
the three m channels. See text for further discussion. 



T'(E). Though the electron hopping between regions 
1 and r of the contact has been neglected in calculating 
the RTM the so calculated channel transmissions approx- 
imate very well those calculated using the FTM so that it 
is very easy to relate the RTM channel transmissions with 
the FTM channel transmission. This shows that the hop- 
ping between the regions 1 and r on both sides of the tip 
atom is almost negligible. Only for the one majority (M) 
channel we see a small deviation near the Fermi energy 
indicating that here 2nd neighbour hopping contributes 
to the transmission of that channel. As the eigenvectors 
of the RTM (see Table HJ reveal, this channel is mainly 
s-typc. Since s-electrons are strongly delocalizcd there 
is a small but finite contribution from second-neighbour 
hopping explaining the deviation between the FTM and 
RTM transmission in that channel. The first minority 
(m) channel is also mainly s-type but now it is hybridized 
with d 3z 2_ r 2 and p z orbitals. The other two m channels 
are degenerate and mainly d xz - and d yz -type strongly hy- 
bridized with p x - and p y -orbitals, respectively. 

As discussed in our previous work^ the five d-type 
transport channels for the m electrons available in the 
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TABLE I: Eigenvectors of the RTM at the Fermi level for the 
contact sketched in Fig. [5] Each column gives the weights of 
the atomic orbitals (AO) given in the left column on the tip 
atom in each eigenchannel shown in Fig. |3] 




FIG. 4: Sketch of Ni nanocontact (27 atoms). As in Fig. El 
the contact consists of two pyramids along the (001) direc- 
tion but now both pyramids share the same atom at the tip. 
The device region (grey circles) consists of 27 Ni atoms and 
the left and right electrodes (empty circles) are modelled by 
Bethe lattices with appropriate tight-binding parameters to 
reproduce Ni bulk DOS. 



perfect Ni chaini^ are easily blocked in a contact with 
a realistic geometry like that in Fig. [2] because the d- 
orbitals are very sensitive to geometry. We have refered 
to this as orbital blocking. It is not so surprising that the 
d x 2_ y 2- and d xy -channel which are very flat bands just 
touching the Fermi level in the perfect chain are easily 
blocked in a realistic contact geometry. These bands rep- 
resent strongly localized electrons which are easily scat- 
tered in geometries with low symmetry. Interestingly, 
even the d3 Z 2_ r 2-channel, which for the perfect chain is 
a very broad band crossing the Fermi level at half band 
width, does not contribute to the conduction as our eigen- 
channel analysis shows. This channel is blocked because 
the d3 Z 2_ r 2-orbital lying along the symmetry axis of the 
contact is not "compatible" with the geometry of the two 
pyramids. On the other hand the d xz - and d yz -channels 
are both open in that geometry because their shape is 
compatible with the pyramid geometry of the contacts. 
This illustrates how the geometry of a contact can effec- 
tively block (or open) channels composed of very direc- 
tional orbitals. Of course, for different geometries we can 
expect different channels to be blocked or opened. 

Obviously, the approximation made in the calculation 
of the RTM becomes worse the bigger the hopping be- 
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(a) majority chan. 1 (b) majority chan. 2,3 (c) minority chan. 1 
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FIG. 5: Transmission functions of open transport chan- 
nels for the Ni nanocontact sketched in Fig. |l]as calculated 
from the FTM T(E) (solid lines) and from the RTM T'(E) 
(dashed lines). The RTM transmissions calculated for the 
cross-section labeled with 0' in Fig. 0]are given by the thin 
solid curves (labeled RTM'). 



tween the regions 1 and r is. For example, in the contact 
geometry shown in Fig. 0] electron hopping from the lay- 
ers immediately connected to the central atom (labeled 
0) is certainly bigger than in the geometry of Fig. 
Indeed, Fig. shows that for almost all channels the 
RTM transmissions differ appreciably from FTM trans- 
missions, making it difficult in some cases to relate them 
to each other. Fortunately, we can judge by exclusion 
which RTM transmission relates to which FTM trans- 



mission since for the other channels at least the RTM 
transmission function mimics the overall behaviour of the 
FTM transmission function. However, for more compli- 
cated situations it might be impossible to match the RTM 
transmission with the FTM transmission for all channels. 
The cure to this problem is obvious: One has to choose a 
bigger cross-section, i.e., add an atomic layer to the cross- 
section so that the hopping between 1 and r becomes small 
again. If we choose, e.g., the cross-section labeled with 0' 
in Fig. 01 (including the atomic layer to the right of the 
central atom) the so calculated RTM transmissions now 
approximate very well the FTM transmissions as can be 
seen in Fig. |5] 

IV. CONCLUSIONS 

In summary, we have shown how to obtain the orbital 
contributions to the eigenchannels at an arbitrary cross- 
section of a nanoscopic conductor. The method has been 
implemented into our ab initio quantum transport pro- 
gram ALACANT and we have illustrated the method by 
exploring the orbital nature of the eigenchannels of a Ni 
nanocontact. The method works very well when the cho- 
sen cross-section is thick enough so that hopping from the 
layers left and right to the cross-section becomes negli- 
gible. Hence in some cases an additional atomic layer 
has to be included to the cross-section we are actually 
interested in. Taking this into account the method has 
no limitations and can be readily applied to ab initio 
transport calculations in all types of nanocontactsi^ and 
molecular junctions. 
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